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\l , We have calculated the single-particle density of states (DOS) for a model of spinfuU Tomonaga- 
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Luttinger liquid with frequency dependent parameter Kc of the charge sector (and Ks = 1 of spin 
sector). Such frequency dependence may originate from interactions with optical phonons. DOS 
exhibits a power-low suppressed asymptotic behaviour near u) = with exponent larger than one, in 
agreement with previous results. For larger frequencies, but still not far away from the origin, DOS 
exhibits a peak, the position of which scales as [i^(0) -|- K~ (0)]aJo, where loo is the characteristic 
phonon frequency. The position of the peak decreases as the phonon's frequency becomes smaller. 
It is interesting that the peak of DOS does not coincide with the characteristic phonon frequency. 

This paper investigates the properties of the single-particle density of states (DOS) at the Fermi surface of ID 
conductors with electron-electron and electron-phonon interactions. Interest in the problem of the electron-phonon ID 
^ , systems has been motivated mainly by many interesting results on organic conductors and superconductors, conducting 
^ ' polymers and also photoemission experiments on the properties of the spectral function near the fermi surface [1], 
, which are known to be influenced strongly by the electron-phonon coupling. The electron-phonon interaction in quasi- 
one dimensional materials has been studied by different authors, but mostly for the case of the noninteracting electrons 
[2-5] . Here the effects of interaction of the charge density with small momentum phonons were treated and the direct 
diagonalization of the Hamiltonian could be done. The combined effects of Coulomb and electron-phonon interactions 
, for phonons with finite frequencies was studied in the case of spinless fermions by [6] and the spinfuU case by the 
' same authors [7]. In these papers the effects of dispersionless phonons with momenta k, 2kp were investigated. In 
O ' both cases the renormalization group approach was used. The effort on those papers was concentrated, however, on 
the effective interaction constants, effective mass due to electron-phonon coupling and renormalized velocities for the 
] charge and spin flucuations. A possible phase diagram (in presence of suitable coupling mechanisms between chains) 
' was represented. 

\ In this paper we consider the effects of electron-phonon, in addition to electron-electron interactions, in the frame- 
' work of Tomonaga-Luttinger liquid model. We focus on large momentum transfer scatterings, of q ~ 2fci?, q Akp, 
ps) ] of electrons by phonons. The electron-electron part of the Hamiltonian [Eq. (2) below] can be written as a sum of 
' two independent (commuting) contributions of charge and spin collective modes. To count for the phonons, in the 
Cp I continuous description, we separate in fast and slow components the deformation field of the lattice sites. We keep 
only the components near 2kp and 4:kp of this field. Physically the 2kF component interacts with electrons in which 
one electron hopes from one Fermi point to the other. Similarily the Akp component interacts with electrons and 
"j^ now two electron can hope from one Fermi point to the other. We integrate out the phonon degrees of freedom and 
^ [ examine the remaining effective electron-electron interactions. The electron-phonon interaction, as observed also in 
I ' the papers cited above, affects the charge as well as the spin sector by renormalizing the coupling constant gs and 
■ O I the Tomonaga-Luttinger liquid parameter of the charge sector Kc- The study here is confined to the case when 

remains negative, i.e. on the case of gapless spin excitations. Our approach is similar to the one of Voit and Schulz 
[6,7], but we differ on the way we treat the retardation effects appearing from the presence of the electron-phonon 
interactions. In [6,7] they use a perturbative expansion to study the retardation effects whereas we investigate their 
effects by studying the effective Luttinger liquid to which the electron-phonon interactions give rise to the frequency 
dependend parameter Kc of the charge sector. We find the imaginary time single particle Matsubara Green's function 
I and numerically produce the single-particle density of states of this unusual Luttinger liquid. 

The model we consider is the one-dimensional (ID) Hubbard model coupled to phonons. The lattice effects can be 
included in this model by making the hopping integral t dependent on the intercite distance: 



o 



o 
o 



X 



t^j ^ t + ^k{u^ ~ Uj) . (1) 



This has been proposed for the first time by Su, Schriefer and Heeger [8] to describe the physics of conducting polymers. 
The Hamiltonian with the above hopping integral has the following form: 
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where Uj is dimcnsionlcss and k has dimensions of energy. The Cj,^ operators are the usual creation and annihilation 
operators for the electrons with spin a in the Wannier orbitals at site j and nj^o- is the number of electrons on site j. 
U is the repulsion of two electrons on the same site. We consider here the case of the incommensurate band filling i.e. 
4fcF 7^ 27r/a, and therefore we do not consider umklapp scattering. 

Substituting the creation and annihilation operators in terms of the new continuous fields i?o- (x) , L„ [x) , and in- 
troducing the scalar U(l) and vector SU(2) currents J =: L\^La. :, J =: Ll^{aa,0/2)L0 :, the first term of (2) can be 
written as a sum Hi = h'^^^^ + Hf^^'^\ Adding to the above part of the Hamiltonian the ?7-interaction term one has 
-H'Hubb = He 



Hg, where 
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j dx {JJ + JJ)+gcj 



dxJJ 
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J dx (JJ + JJ) + 5s y 



da;JJ 



(3) 



with Fermi velocities Vc = vp — Ua/i^iirvp)) and = Vpi}^ — ^Ua/ {^iirvp)) and coupling constants Qc = Ua/2, 
Qs = —2Ua, in the charge and spin sectors respectively. The coupling constant in the spin sector is negative which 
results in the spin excitations being gapless (the ggJJ term is marginally irrelevant in the case Qs is negative). Using 
Abelian bosonization one can restore to He given by (3) the canonical Gaussian form: 



Vc 
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dx 
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(4) 



whith Vc and depending on the coupling constant gc in the following way: 77c = Wc(l — 9c /''^'^'^c Y^^ — [(1 ~ 

(7c/7i'Sc)/(l+.9c/7''Cc)]^^^- We will see that the electron phonon interactions lead particularly to strong renormalization 
of the parameters of the charge sector. 

The cifFect of electron-phonon interactions on the dynamics of the electronic system can be easily calculated since 
the pliononic action is quadratic in the displacement field whereas the clcctron-plionon interaction is a linear function 
of the lattice displacement (multiplied by the electronic staggered charge density). Integrating the displacement fields, 
in the integrals for the partition function, one gets in addition to (3), an effective interaction between the staggered 
components of the electronic charge densities (see for instance [9]) 



Si^t = -2TrVcXi fdrdT'dxl ^ f d^xdr' p*i2lkF,x,T') (5(r' - r) + -^^^^e"'^'!^'-^! 



p{2lkF, X, t) 



(5) 



where 



p{2kF,x) = ^^Rl{x)La{x) = Ocexp |iV27r<i?c(a^) +2ikpX | cos{A/27r$s(a;)} , 



(6) 



p(4A;f, x) = ttu exp i^iVSn^dx) + 4ikpX j , (7) 

are the 2kF and 4fcF-components of the staggered charge density end wi = w{2lkF) for I = 1,2 are the phonon 
frequencies and A; = [Ksm{lkFa)]^ /{TTVcPiUifa'^), where pi is the effective mass of the l-th mode. It is of course known 
[and can be seen also in the above formula (5)] that the electron-electron interaction mediated through phonons is 
constituted of two kinds of interactions. First there is an instantaneus attraction effect between the electronic staggered 
charge densities. By observing the formulae for the charge densities Eqs. (6-7) one can see that only the 2fcF component 
gives a contribution to the instantaneous component of the interaction. The 4A;f component of the staggered charge 
density have the particular exponential form of complex argument and does not give any contribution. Quantitatively 
the effect of the instanteneous attraction between the 2A;f components of the staggered charge density can be evaluated 
using the identity [10]. This shows that this kind of interaction renormalizes the LL parameters of both the charge 
and the spin channels. In the latter the additional constant will be with positive sign and will tend to increase to 
a value which we will assimie to be still negative. A positive Qs will lead to a gap in the spin sector, something we are 
not undertaking to examine here. 

The contribution from the nonlocal term, or the retarded attraction between the staggered charge densities, will be 
discussed a bit longer in the following since they are the factor that leads to the frequency dependence of the Luttinger 
liquid parameter Kc in the charge sector. Let us see first the effect of the interaction between the 4A;f components of 
the staggered charge density: 
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J dx dr' : p{4kF,x,T) : (^^^re""'"'^"^'') : p(-4fcF, a;, r') : 
= - fdxdr' : drP{4kF,x,T) :: dr'p{-^kF,x,T') : Z^^g-'^^l^-^'l 
= -Aw^cal —drMr) j dr'e-'^^l^'l exp{iv^[$e(r) - $c(r - r')]}5^$e(r - r') 

W2 J 



If we note with 



1^2 J 



(8) 



in the following we make the mean field approximation by representing the Akp contribution in w-space in the following 
form 



j dT[a^$c(T)] [drFikAT-)]= j dujJ^<^c{-oj)FikA'^) = j dujoJ^^c{-oj)hkF{<^)^c{'^ 



(9) 



The 2kF components give for f2kF two terms from which, within this approximation, only one gives a finite contribution 
in the dynamics. The first one being 



hakAr') = TT^t^ea^ -e-^l-'l (ei^*^(-)e-i^*=(---') cos [v^$.(r) 



COS 



27r$,(r-r') ) 



is similar to fikp ■ The additional one is 



/2,2fe,(r') = 7r2t;,a2 ^e-il-'l (ei^*=(-)e-i^*=(^-^') ) 
multiplying 

dr COS ^/2TT^a (t) dr COS 

In w-space the term (10) will give a contribution similar to (9), namely 

dT[dr^,{T)][drF,,2kAr)]= j duj u'' ^ ,{-u)h^2kA^)^ M 



27r$s(T-T') 



(10) 



(11) 



(12) 



(13) 



whereas the term (11) can be shown to be a total time dcrvativc and therefore wc do not add it as a contribution. 
In the contributing terms f^kp and /2,2fcF we have made the mean field approximations 



and 



exp {iV8^ [$e(r) - $c(r - r')]} « (exp {iVS^ [$e(r) - $e(T - r')]}) = (^^) 



exp |i\/27r [$c(t) - $c(t - r')]| cos V2tt<^c{t) cos ^/2tt^c{t — t) 



(14) 



(exp{i^[$e(T) - $c(r - r')]})(cos 
where d4kp = and d2kp = = 1/2 and 



2d2 



^$e(T) COS y2^$e(T " t') ) ~ ( ^ ) " ( ^ 



2d, 



2ro 



4is:o - 1 



(15) 



(16) 



To is a low cut-off of the order of Cp^. The integral (16) is always positive since in the Hubbard model takes values 
decreasing from 1 to 1/2 as u=U/t increases from to oo (see for instance Frahm and Korepin [11]). 



In cj-spacc the effect of tlie nonlocal term of (5) is that the charge sector can be described by the modified Gaussian 
action with the kinetic energy nonlocal in time; 



where the function /(w) takes values between 

/(w » 0) := 1 + v,Kl J dr' [fi,2kAr') + hkAr')] = 1 + 2Trv^,K°To ja^^ + 4al 



102^0 - 1) 



(17) 



(18) 



and f{uj 3> 0) :~ 1. In Eq. (18) we introduced the w-independent constant k. 

It is interesting to explore the consequences of frequency dependence of /(w). Its form, within our approximation, 
is given by 



fico) = 1 + 



(19) 



where /(O) = k'^ (the parameter k takes large values) and /(w » 0) = 1. wq is considered to be the phonon's frequency 
(the phonons are considered as without dispersion and we can approximate uj2kF ~ ^ikp = ^o)- Large values of k 
notify strong renormalization, and this is stronger at smaller frequencies (but not smaller than the critical one set from 
the spin sectors requirement of no gap opening). The value of the constant k is expected to be big, and experimental 
observations in materials like Ko.sMoO.i [3] give a value of several hundred. 
In time representation the action on the charge sector has the following form 



da;dr 



Vr 



+ drMr) J dr' [fi,2kAr') + UkAr')] dr^r - t')| • 



(20) 



The effective retarded attraction between the staggered components of the charge density affects only the charge 
sector. In the spin sector remains 1, [preserving the SU(2) symmetry]. 
For the special limit a; = in f{u)) of (17), it takes the following form 



-k^ {dr^.f + V, (5,$c)' 



[/('^ = 0)] , 



which in the canonical form (4) is described by the new renormalized parameters 



The limit w in the charge sector is described by the action 

1 
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dxdr 



Vc 



[/(^»0)] . 



(21) 



(22) 



(23) 



We consider the case of being near 1 and neglect the renormalization coming form the local term, which is of the 

same magnitude as the renormalization of the spin sector and which we considered to be small. In this case Kc = l/k. 
The single-particle Green's function is expressed as a correlator of the chiral exponents of the charge and spin fields: 



• exp 



iV27r(/)c(0,0) 



G{t,x) = (exp iv27r(^c(''', x) exp i V 27r(/)s (r, x) 

= (r — a;/us)~"^/^(exp i\/27r(/)c(T, x) exp i-\/27r0c(0, 0) ) 
The model (17) leads to the Matsubara Green's function 



exp 



iV27r<^s(0,0) ) 



G{x = 0;t) = - exp < -- 



da; 



1/4 



sin^ {lootx/2) 



(24) 



(25) 
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This Green's function for large r behaves as 

G(r) =^ p^^A.)^-S^.^^ ; ^c = i(^-V^)' ; [/(u. = 0)] , (26) 

which corresponds to /(cj = 0), Eq. (21). This is so because the exponential of the Green's function Eq. (26) for large 
r behaves as — [(v^— l/-\/fc)^/4] Inr. 
For small r, it bahaves as 

G(t)~^ =^ p^M^'^-^Konst , (r^O) ; [/(a; » 0)] , (27) 

corresponding to the case with /(w) with large w, Eq. (23), [if G{t) l/r were in all r space then G(ia;„) = 

~2i J^°° dr sm{ijJnT) /t = — 2iSi(+oo) = — itt and p^^^^{lo) = {—1/tt)QilrG^ = 1]. 

To extract DOS on the whole frequency space one has to Fourier transform this function to get G(iWn) and then 
take the analytic continuation icOn ui + iO. The values of p^^^ (u) between the two asymptotic limits, w — > 0+ and 
a; ^ 0, arc shown in Fig 1. They arc obtained numerically by analytic continuation of G{iLUri) from the imaginary axis 
to just above the real axis. A remarkable feature of DOS is the peak at w « K~^u!o/2 (see Fig. 2). 

Since DOS is invariant under K — > this empirical formula can be generalized as 

ojpeek = ^{K-^ + K)uo . (28) 

It is interesting that the peak always occurs at frequencies larger than the characteristic phonon frequency wq- At 
small (large) values of K this discrepancy can be quite substantial. For example, a peak in DOS has been observed in 
K0.3M0O3 at w w 300 meV. The behavior of DOS at smaller frequencies is almost linear in w which suggests K k. 0.15. 
Then Eq. (28) gives a reasonable estimate for the phonon frequency: ~ 90 meV. 

For smaller the Luttinger parameter becomes small, the exponent 9c becomes large and the peak approaches 
the origin. This is also observed in the numerical plots, and is in accord with (28). 

The Fig. 1 clearly show also the crossover from the Luttingcr-liquid type behavior with singular dp/duj &t K > 0.17 
to the power-low suppressed with exponent larger than 1 behavior &t K < 0.17. 

To summarize, we have investigated in this paper the single particle density of states at the fermi level of ID 
conductors at inc;oninicnsuratc fillings with clcctron-phonon interactions. The problem is studied using the Tomonaga- 
Luttinger liquid approch. The phonon degrees of freedom were integrated out and the remaining effective electron- 
electron interactions were studied. The electron-phonon interactions with large momentum transfer of g w 2kF and 
q w Akp result in retardation effects which lead to the unusual form of LL model, namely with parameters of the 
charge sector depending on the frequency. The LL constant K is also a function of wq, the phonon frequency. At small 
frequencies Kc and Vc are strongly renormalized by the interactions. We found the imaginary time Matsubara Green's 
function for this model and built the single particle density of states. Dos exhibits a peak which is proportional with 
the phonon frequency, but nevertheless they do not coincide. 

E.P. is grateful to Steve Allen and to Dave Allen for helpful discussions and to A. M. Tsvelik for suggesting the 
problem and discussions. 
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